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SUMMARY 


This report deals with methods of measuring the probability distri- 
butions and mean values of random functions as encountered in turbulence 
research. Applications to the measurement of probability distributions 
of the axial velocity fluctuation u(t) and its derivative du/dt in 
isotropic turbulence are shown. The assumption of independent proba- 
bilities of u(t) and du/dt, which has been used as an approximation 
in the application of zero counts to the measurement of the microscale 
of turbulence A, is investigated. The results indicate that the assump- 
tion is satisfied within a few percent and that there is, so far, no 
evidence that the systematic difference between A measured from zero 
counts and A measured independently can be traced entirely to the sta- 
tistical dependence of u and du/dt. 

The chronological development of apparatus is described, concluding 
with the present 10-channel statistical analyzer based upon a system of 
pulse amplitude modulation followed by an amplitude discriminator and a 
counter. A discussion of the relative merits of various systems is 
included to indicate the reasons for this choice. 


INTRODUCTION 


Let l(x,y,z,t) represent a quantity such as a velocity component 
or a pressure in turbulent flow. Since turbulence is an essentially 
statistical phenomenon, l(x,y,z,t) will not be representable in the 
same way as an ordinary function but will be defined only by certain 
probability distributions and mean values. Experimental studies of tur- 
bulent flow are thus primarily concerned with the measurement of mean 
values. There are several ways of defining the mean values of l(x,y,z,t) . 
Often the most convenient, theoretically, is to define I as an ensemble 
average; that is, one considers N similarly prepared systems, say wind 
tunnels with the same grid setup, and measures l(x,y,z,t) simultane- 
ously at corresponding points of these N wind tunnels. The N results 
can then be evaluated statistically and yield the probability distribu- 
tion and the desired mean value. 
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Experimentally, and also often in theory, one uses instead the time 
average. One observes l(x,y,z,t) for a sufficiently long time in one 
system and defines 


provided this is possible. If there exists a time -independent mean 
value, the integral will reach a constant value for a sufficiently large 
value of T. Clearly one can also use a space average, defined by 


where V is a certain volume in space. Finally, one can combine the 
latter two and average over both space and time. 

Most of the recent experimental work in turbulent flow is con- 
cerned with flows which are steady in the mean, that is, for which the 
time average exists and for which the time and ensemble averages can 
be rigorously interchanged. Investigations of isotropic turbulence 
behind grids and of wakes, jets, channel flows, and so forth belong to 
this group. Isotropic turbulence set up in a box, on the other hand, 
is clearly a problem for which a time average has no meaning: here, 

evidently, the ensemble average or space average is wanted. 

The mean value of a power of I, for example I 2 or I^, can be 
measured in essentially three ways. First one can use a device which 
squares or cubes the instantaneous value of l(x,y,z,t) and then pro- 
ceed to take the desired mean value. Secondly one may use a device 
which produces the mean of a certain definite power directly, such as 
a thermocouple meter for measuring the time mean square of a current. 
Finally, one may obtain the probability distribution or probability 
density p(|), where p(|) d£ denotes the probability of finding I 
between £ and £ + d£, and then obtain all the mean powers of I as 
the moments of p(|). Whether a time or ensemble mean value results 
depends upon the determination of p(£). If p(|) d£ denotes the 
fraction of time the function l(x,y,z,t) spends between £ and 
£ + d£ divided by the total time T, a time average is taken. If 
p(S) d| is defined as the number n of systems having values of I 
between | and | + d£ divided by the total number of systems, then 
an ensemble average is taken. 




V 


NACA TN 5037 


3 


The methods described below are primarily intended for the measure- 
ment of time averages and are based on this last form of the averaging 
procedure. The original apparatus handled p(|) directly, but it was 
subsequently found more convenient to work with the integrated proba- 
bility density, or distribution function, F(|), which is defined as 
the probability of finding l(x,y,z,t) with a value greater than £ 
and is related to the probability density as follows: 


p 00 

F(S) = / p(£) d£ (1) 

J £ 

Clearly F(-a>) = 1 and F(°°) =0. It follows that 


dF = -p(g) d| (2) 


Although time averages are dealt with here, it is also possible to 
use the same procedures for ensemble averages. 

The apparatus was built with the idea in mind of reducing the prob- 
lem of finding the probability density, or the distribution function, to 
one of pulse counting. The advantages of such a procedure are that 
counting methods permit averaging over almost unlimited intervals of 
time, and pulse-counting circuits are in general very stable and satis- 
factory. This approach grew out of the investigation of the zeros of 
fluctuating-velocity components (ref. l) , reported earlier, which 
demonstrated the convenience of pulse-counting methods . 

This work was carried out under the sponsorship and with the finan- 
cial assistance of the National Advisory Committee for Aeronautics. 

The authors wish to express their appreciation to Mr. George Skinner 
for his generous help in the preparation of the manuscript. The original 
gate apparatus was designed and built by Mr. C. Thiele and Dr. H. Wright, 
who furnished very useful advice during the early stages of the work. 

Part of the work reported here was done as early as 19^9 • The present 
report summarizes both the earlier work and recent developments (ref. 2) . 


SYMBOLS 


t time 

amplitude coordinates 
u(t) axial velocity fluctuation 
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V(t) any fluctuating input voltage with zero mean 

~K microscale of turbulence 

p(|) probability density 

p(|,T]) joint probability density 

F(|) distribution function 

Dimensionless parameters: 

C = 



S = 


= Skewness 


f = 


(^) 2 


= Flatness 


For a Gaussian distribution: 


c = # 

S = 0 

f = 3 

PRINCIPLES OF APPARATUS 


In all cases the apparatus handles a fluctuating voltage V(t) 
(figs. l(a) and 2(a)) such as the output of a hot-wire and amplifier 
combination. The term V(t) may represent a velocity fluctuation 
measured at a certain point, its derivative, or so forth. The original 
experimental apparatus was set up to determine the statistical prop- 
erties of V(t) by measuring its probability density step by step, 
using a gate circuit essentially similar to Townsend's "statistical 
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analyzer" (ref. 3 ) which passes a signal only if the amplitude of V(t) 
lies within a narrow range of voltages. Specifically this circuit pro- 
duces a signal of constant amplitude during the time that V(t) lies 
between £ and £ + d£ , where £ is the voltage setting of the gate . 
Hence the operation performed by the gate circuit yields a signal like 
the one shown in figure l(b), consisting of irregularly spaced square 
waves of uniform height. The total area of these square waves, or the 
s\om of all their widths, measured over a time interval T, thus gives 
the time T^ which V(t) spends between £ and | + d| during T. 

The probability density p(|) d| is defined by stating that p(£) d| 
denotes the fraction of T which V(t) spends between £ and £ + d| ; 
that is, 

p(|) d| = ^ (3) 

T 


To measure the sum of the widths of the square segments of fig- 
ure l(b), a counting method is used as follows: A pulse generator pro- 

duces equidistant pulses at a rate of, say, M per second (fig. 1(c)). 

A coincidence circuit mixes the signals in figures 1(b) and 1(c), giving 
the result shown in figure l(d) . The signal in figure l(d), which con- 
sists of a set of "chopped" square waves, is fed into a counting circuit. 
Let m be the number of counts obtained per second, then evidently 


pU) « » S 

M 


( 4 ) 


Varying £ then yields p(|) completely. From p(|) the averages 
are obtained as the moments; for example. 




l 3 p(l) d| 


p(t) d| 
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The time over which a count is taken is evidently limited only by the 
stability of the whole setup. Especially important for large values 
of £, for which p(|) is very small indeed, is the fact that the time 
of measurement can be increased and the accuracy very much improved. 

This first method lends itself well to direct measurements of mean 
values by making use of Simpson's rule in determining how long to let 
the apparatus run at each value of | . This technique is described in 
the section "Obtaining Mean Values." 

To reduce the measuring time involved, a number of gate circuits 
can be combined into a multichannel analyzer, so that p(|) can be 
completely determined in one short run. However, when a multichannel 
arrangement is considered, it becomes evident that one should change to 
a system which measures F(|) rather than p(|). The reasons for this 
will be discussed in the sections "Details of Apparatus" and "Errors 
in Apparatus . " 

The principles of this multichannel form of the apparatus are 
illustrated in figure 2. Since the counting circuits require their 
information in the form of narrow pulses, the incoming signal V(t) 
is first made to amplitude modulate a continuous train of pulses 
(fig. 2(b)) in such a way as to result in a pattern (fig. 2(c)) which 
lies completely above the zero voltage line. A discriminator, or rec- 
tifier, which is biased to a voltage value of | is then used to cut 
off any pulses whose height does not exceed | (fig. 2(d)). Suffi- 
cient amplification of those pulses which pass the discriminator is 
provided to insure that each will register in the counting circuit 
which follows. In the multichannel setup, 10 discriminators are used, 
each set at a different | level. These levels are adjustable to per- 
mit the use of all 10 channels for any portion of the incoming signal 
to obtain, for example, 10 points on the distribution function for 
positive values of V(t) and 10 points for negative values of V(t) 
by reversing the input signal. 

It will be observed that in passing from V(t) (fig. 2(a)) to the 
train of modulated pulses (fig. 2(c)) the zero level of the signal has 
been lost. This must be redetermined later using the knowledge that, 
by definition, v(t) is a fluctuating voltage whose mean value is zero, 
giving a value £ 0 which corresponds to V(t) = 0. One can then say 


6 -5 - Cq 


( 5 ) 
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If, again, the pulse generator supplies a train of pulses at the 
rate of, say, M per second (fig. 2(b)) and m is the number of counts 
obtained per second, then evidently 


F (d - = 


( 6 ) 


OBTAINING MEAN VALUES 
Direct Measurement 


With the original gate circuit for measuring p(|), it was possible 
to measure mean values directly, by the application of Simpson's rule; 
that is, the integration giving the appropriate moment of p(| ) could 
be carried out with the apparatus itself. 

Assume that the pulse generator furnishes M pulses per second and 
that for a given value of | the count is taken over t^ seconds 

leading to a total of m» pulses counted. Then: 


p(6) d| - -i- 
Mt. 

or 

t|p(l) dl = Hi (7) 

s M 

Assume now that the zero moment that is> 





p(s) 


is to be measured. Using Singpson's rule the integral can be written 
as the sum 



p(S) 


i|p(0) + 4p(l) + 2p(2) + 



(8) 
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±f -the interval | n - I n-1 -*- s e Q. ua b b bor convenience. Similarly 

the other part of the integral may be written 



p(i) d£ = l[p(0) 


+ 4p(-l) + 2p(-2) + 


•] 


(9) 


Since p(|) goes to zero for large values of | £| it is irrele- 
vant how the last point is chosen, provided the sum is extended to large 
enough values of | £ | . 

Consider equation (7) and choose t| in such a way that t Q = t, 
t +1 = At, t +2 = 2 t, and so forth. Then clearly 



”1 


= Mt p(0) + 4p(l) + 2p(2) + 


•] 


that is, 


and similarly 


Z 

1=0,1,. 


me = 3 Mt / p(t) d£ 
J 0 


z 

t = 0 ,- 1 ,... 



p(i) dl 


(10) 


( 11 ) 


Thus if one sets the device successively to £ - 0, 1, 2, 3, • • • , 
counts time r, 4 t, 2t, . . . seconds, and leaves the numbers in 
the register, the apparatus performs the summation and by virtue of 
equation (10) the integration for positive values of I . By reversing 
the input the integration for negative values of | is obtained. 

The procedure can be applied in principle to any moment inte- 
gral Oh, since 

| n p(l) d£ n = 1, 2, 3> Z • 
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Now t| is chosen as follows: 


t Q = ° 


Thus 


t|p(|) = 


= 4t( ± l) n 

t+2 = 2r(±2) n 

" 21 * “pU) 

n = 2,4 I 

4t|Vi) 

n = 1,3 J 

2t(-1)Vp(|) 

( 

-=fr 

CVJ 

II 

g 

M-i) n s n pU) 

> 

n = 1,3 


I >0 


£ <0 


( 12 ) 


and thus 


= 3Mt / | n p(|) d| 


(13) 


l>0 


and 


^ = (_i) n 3Mr j tvo di 

?<0 

A procedure like this is sometimes convenient. Since the functions 
are smooth, not very many points are required to approximate the inte- 
grals very closely. It is also evident that this procedure could well 
be performed automatically if the additional apparatus were warranted. 


Short Cuts 


For some applications it is possible to shorten the procedure very 
much by using a known form of p(|). For example, assume that V(t) 
represents a velocity fluctuation u(t) in isotropic turbulence. In 
this case one knows that p(|) is very nearly a Gaussian distribution. 
Thus: 



p(l) = p(o) e 
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and hence 


p(£ x ) = p(0)e 


1 5 1 




Thus 


u 


2 




2 log 


p(0) 

P(^l) 



( 1 ^) 


and thus is obtained as the ratio of two measurements. Similarly 

can be obtained. 


Computation of Mean Values Using Multichannel Analyzer 
Consider the mean value of the nth power of the input V(t) 


V n = 


SV6) d| 


p(S) 


Now 


p(6) d| = 1 


so that 


V n = 


l“pU) d| 
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and by equations (l) and ( 2 ) 


V n 



i n dF 


Because of the factor £ n , the integrand becomes infinite at the limits, 
making it unsuitable for numerical evaluation. This difficulty can be 
overcome by transforming to an integral with respect to the |-axis. 
Splitting the integral into two parts about | =0: 


V n = [ £ n dF + f £ n dF 
J 0 J a 


pa nO 

/ £ n d F - / £ n d(l - F) 
'0 J l-a 


where F(0) = a. Integrating by parts: 





n| n_1 (l - F) d£ 


or 



8 

II 

WLf 

r 

L 

p°° 

P-oo 

U n (l-F) 

r ^ + 

I n_1 (l - F) d£ 

L i=0 

Lyo u o 


The first two terms vanish at the limits if F(n) approaches zero 
exponentially or faster as | approaches minus infinity. This require- 
ment is always satisfied if the mean value is finite. Therefore 


V 


n 


n 



| n_1 F d£ + 



| n_1 (l - F) d£ 


(15) 
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The Integrals are evaluated numerically by Simpson's rule. (See fig. 3 .) 
Since the main concern in turbulence work is distributions approaching 
a Gaussian form, the fourth moment of a Gaussian distribution was com- 
puted by Simpson's rule, using nine points. The result deviated from 
the theoretical value by less than 1 percent. 

In practice, the normalization procedure 


poo 

/ p(!) dg = 1 


involves dividing the number of counts m^ in any register by the total 

number of pulses generated by the pulse generator during the counting 
time Mq so that: 


V n 



m e 

I 11 - 1 _L d| + 


Mq 





(16) 


One can conveniently run the generator for 1 minute with a pulse rate 
of 10 ^ per minute, so that 


*0 


m^_ 

io5 


Finally, it may be called to mind once again that the exact posi- 
tion of |=0 is not known beforehand. Using the fact that the first 
moment of the distribution function is equal to zero (that is, V = 0) , 
the position of this axis may be determined. This procedure has proved 
more satisfactory than any in which one tries to set the initial "zero" 
pulse amplitude to some known arbitrary value with the accuracy required 
for confuting mean odd powers. 


DETAILS OF APPARATUS 
Original Gate Circuit for p(|) 


The circuit diagrams of the elements are shown in figures L and 5 . 
The block diagram is given in figure 6 , which should be self-explanatory 
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The mixer will he seen to he a twin triode normally running saturated 
through a common anode load resistance. Strong negative grid signals 
will cut off the triodes, hut a large anode signal will appear only when 
hoth halves of the tuhe are cut off simultaneously. The small effect 
of one-half being cut off alone is omitted in the block diagram. It 
should he pointed out that one-half of the first mixer receives a posi- 
tive square wave which causes grid current to flow, thereby building up 
a steady negative bias nearly equal to the peak voltage of that square 
wave, so that it considers itself to be receiving a negative input 
during the periods when this positive square wave is absent. As the 
5 + d£ line approaches the peak of the input signal, this positive 
square wave becomes of very short duration and it is difficult for the 
mixer to hold its bias. In operation the steps shown in figure 6 must 
be studied on an oscilloscope and levels adjusted, so that the sundry 
effects, such as those mentioned above, do not interfere with the 
counting of the desired pulses. 

The signal is first amplified to about a 100-volt peak, to enable 
the gap to be of the order of a few volts and still be small compared 
with the peak amplitude. The Helipot A (fig. 4) adjusts the triggering 
point of the gate, that is, 5; the two resistances B and C adjust the 
operating points of the two sides of the gate independently and hence 
set the gap width, that is, dg . 

Adjustments are made by studying the output of the entire device 
on an oscilloscope, while a sine-wave signal is being fed in. The posi- 
tion of A for 5=0, the width of the gap dg in volts, and the ampli- 
fication factor giving 5 or the setting of A in volts can be determined 
in this fashion. 

Pulses can be generated in any convenient manner, for example, by 
differentiating a square wave and suppressing the positive pulses pro- 
duced, but for many applications it is quite good enough simply to use 
a sine wave. The frequency was ordinarily chosen to be 20 kilocycles 
per second, but different frequencies have been used. 

The resulting pulse sequences are counted on a binary scalar counter 
feeding a mechanical register having a maximum rate of 60 per minute and 
capable of registering up to 10^. The binary system is standard practice 
and uses diode coupling. The mechanical counter (a Cyclotron Specialities 
product) can be switched in to count every 2^ to 2^ pulse by tapping onto 
the output of the fifth, sixth, . . ., ninth stages of the binary scalar 
counter. This switching is useful in obtaining more counts on the register 
at low counting rates (in preference to reading the neon light of each 
stage and adding to the register reading) and may also be employed when 
integrating by Simpson T s rule with the apparatus. 
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Multichannel Analyzer 

In the section ’’Principles of Apparatus” the basic principles were 
outlined and these are shown pictorially in figure 2. The incoming 
signal V(t) (fig. 2(a)) is first added to a train of pulses, having 
a rise time of about 2 microseconds, obtained by differentiation of a 
square wave of frequency up to 10 kilocycles per second. Figure 7 
shows an oscillogram of a single pulse, while figure 8 shows the sum 
of V(t) and the pulses, when V(t) is a sine wave. This step is 
accomplished by means of a differential amplifier (ref. 4), which 
handles an input of about 1 volt peak to peak across l/2 megohm, at 
frequencies ranging from l/2 to 50,000 cycles per second. By means of 
rectification and direct-current restoration, the lower part of the 
signal (fig. 8) is removed and the line A-A is made the zero level, 
leading to the result shown in figure 9, which has a maximum level of 
about 90 volts . The circuit diagram of the pulse amplitude modulator 
is shown in figure 10. 

The modulated pulses are then fed into 10 amplitude discriminator 
channels each consisting of a biased diode (fig. ll) which passes that 
portion of any pulse which exceeds the bias voltage. For example, a 
51-volt pulse will pass a discriminator biased to 50 volts as a pulse 
of 1-volt amplitude. The amplifier which follows each discriminator 
is capable of causing a 0.05-volt pulse to be counted, which corresponds 
to a sensitivity of better than 1 percent of the bias difference between 
adjacent channels. At present the bias levels are set by dry cells, but 
it is intended to replace these by a regulated power supply and a chain 
of precision resistors. To obtain intermediate bias settings (for 
example, if it is desired to double the number of points read, by 
making two runs) the diode bias ground tap may be shifted. 

To minimize the effect of differences in the characteristics of 
the diodes they are heated by 3-9 volts instead of the usual 6 volts. 

This also reduces the quiescent plate currents. 

Each channel is provided with from one to four commercial electronic 
decade scalars followed by a five-place impulse-type mechanical register. 
The decades operate on the ”f lip-flop” principle (ref. 5) and indicate 
their count by means of neon lights. The mechanical register is coupled 
to the last decade by means of a cathode follower and power amplifier. 

A maximum of 10^ counts may be stored in each channel. 

The analyzer may be used to study all, or only part, of an incoming 
signal according to where the modulator output is set relative to the 
discriminator voltages. When dealing with a complete distribution 
function, the input is disconnected and the modulator output (uniform 
pulses) is adjusted so that the middle channel just starts to count. 
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This, as has been explained, is not used as the zero level, which is 
determined later by computation. To count only one-half of the signal 
at a time, the same procedure is followed, except that the output is 
adjusted so that the lowest channel just starts to count. By reversing 
the input signal, the other half may be analyzed, giving a total of 
20 points on the distribution curve. 

A run is made by switching the modulator output to the amplitude 
discriminators for a measured time. For a 1-minute run the timing 
error is less than l/2 percent. 


ERRORS IN APPARATUS 
Stability 


In the original gate-circuit setup, stability was exceedingly 
important because of the measuring times involved (up to 2 hours in 
some cases) . Long "warm-up" periods were essential, particularly to 
insure the stability of the | = 0 level, as this is of paramount 
importance in determining means from the odd moments of the probability- 
density curve since this involves a small difference between large num- 
bers and is extremely sensitive to errors in the | = 0 setting. For 

example, the most important skewness factor in turbulent flow (du/dt)^ 
is obtained by measuring the distribution of du/dt, p(q) say, and 
forming the integral 



dTj 



t] 5 p(tj) dt) - 



t]3 p (t]) di] 


either by computation or in the apparatus. The real difficulty arises 
because the relative position of the zero level of the incoming signal 
and the zero level of the gate or multichannel discriminators may change 
with time. In using the 10-channel analyzer for measuring F(g), it 
was estimated that, in the case of the above skewness factor, a 3-P® 1 ' ce nt 
error would arise if a constant error in the position of the | = 0 level 
of l/2 percent of the interval between adjacent channels existed. 

Generally one can achieve sufficient stability, particularly for the 
short times required for a run with the multichannel apparatus, and by 
computing the £ = 0 axis from the first moment of the probability dis- 
tribution or by equating areas of each half of the distribution function 
(see appendix) the error can be minimized, in the case of the multi- 
channel setup, to 0.1 percent of the interval between channels. 
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Hysteresis 


Any circuit designed to trigger at a certain voltage will exhibit a 
certain amount of hysteresis. Some hysteresis is necessary for stability. 
In effect, the circuit triggers at one voltage when the input is increasing 
and at a slightly lower voltage when the input is decreasing. The effect 
of this imperfection can be reduced by increasing the voltage of the 
input signal. However, too large an input signal results in a very short 
time of transit through the gap, in the case of the original gate circuit, 
and the final output square waves become distorted, making it impossible 
to count properly. In the case of the multichannel analyzer, the use of 
biased diodes as discriminators has reduced hysteresis to a negligible 
amount. 


The finite width of the gap does introduce a correction, which, 
however, can be estimated. If the gate width is equal to e, say, then 
what is really measured is the quantity 


Gap Width in Original Gate Circuit 


1 

e 


l 


l+e 


^ ^ ^meas 




p(t') 


(17) 


and not p(|) d£ =p(0 € - For a small e, however, one evidently has 


I. 


S+€ 


p(l') d|' = p(l)e + p'(|) ^ + . 


p(l)efl + 



I 


and for a known distribution the correction can be evaluated. For 
example, if p(l) is nearly Gaussian 



(18) 
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Ordinarily measurements have to be taken up to values of | of 

the order of Hence the maximum error near the outer edges of 

the distribution is of the order 

meas ^ 3 € 

p(h ~ ' 2 {p 

This error may be appreciable if fairly large gate widths are used. 
However, the errors in the dimensionless ratios of moments of p(|), 
namely, the skewness and flatness factors, are much smaller and for the 
most part negligible. For example, the errors in the flatness factors 
can be evaluated from 



p(S) 


meas 


p(i)[i-4 

21 2 


The flatness factor f is 

f iVt) at 


f = 


sSU) di 


2 


p(l) d| = 


(19) 


U 2 ) 1 


( 20 ) 


The factor f obtained from the measured distribution can be written 
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In the correction term a Gaussian distribution can be used and since 


f |V^ 2 « 

] 

r( n + lN | 

‘0 


\ 2 J 


= P- 


r e^-p 2 e 2 d6 

r (f) 

'o 



( 22 ) 


it follows that 


f 




+ 



Hence even if the gate width e is as large as 20 percent of 
error in f is small. 



the 


Values for Multichannel Analyzer 


As has been mentioned, the amplitude discrimination is reliable to 
0.05 volt while the maximum pulse amplitude is about 90 volts. The 
amplifiers driving the decade counters can be trimmed for gain, so that 
when all are tied to the same level their readings will not deviate more 
than l/2 percent from each other. Without trimming they are already 
good to 1 percent. 


Using dry cells to set the discriminator levels, a spacing of 
yi. volts between channels can be maintained uniform to 1~ percent accu- 
racy with negligible drift over a period of a few months. 

The zero level, as was explained earlier, can be computed to 
0 . 1 percent of the interval between channels . 

Over-all repeatability has been checked for a sine-wave input and 
found to be ±0.1 percent of the reading at each level. Measurements 
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of du/dt in a turbulent air stream gave the following scatter of the 
mean values computed from consecutive 1-minute runs: 


V I 2 +1 percent 


x: +2 percent 


V 5 ~ ±3 percent 


~ ±4 percent 


Main Improvement Introduced by 10-Channel Setup 


It should now be clear that, whereas the gate circuit draws its 
information from selected narrow bands of the incoming signal, the 
10-channel device utilizes the whole signal. This can, of course, be 
done with the gate circuit by using only one trigger circuit and suit- 
ably rearranging the block diagram of figure 6, but, because of the 
long time involved in measuring, the stability requirements become pro- 
hibitive. Operating 10 channels as a p(|) device with large intervals 
defeats the advantage of using p(£), namely that moments can be com- 
puted in the apparatus by suitable timing (see section "Obtaining Mean 
Values"), since the correction outlined above in "Gap Width in Original 
Gate Circuit" has to be applied and therefore leads to computation work. 
It will be realized, too, that measuring F(|) instead of p(|) leads 
to certain simplifications in circuitry since each channel has to dis- 
tinguish only whether each pulse does or does not rise above its bias 
level. These were the considerations which led to the adoption of the 
present arrangement. 

Other schemes will be found in references 6 to 9, but the present 
arrangement is designed for relatively high counting rates. 


The multichannel analyzer was checked out on a sine wave (190 cycles 
per second) . The results for one-half of the signal are shown in fig- 
ure 12, together with the theoretical result, and show close agreement. 


Figure 13 shows the results of measurements of p(|) with the gate 
circuit, for the axial velocity distribution in isotropic turbulence 
behind a grid. The Gaussian distribution is shown for comparison. 


SAMPLE APPLICATIONS 


Sine Wave 


Isotropic Turbulence 
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Agreement is good and confirms observations by others, especially 
Townsend. Skewness and flatness factors are shown in the table in the 
next section. 

Figure 14 shows F(|) for the same quantity under practically the 
same conditions, as measured on the 10-channel analyzer. 


Isotropic Turbulence, du/dt and d 2 u/dt 2 

Figure 15 presents the results of measurements of the distribu- 
tion p( T]) of the velocity derivative du/dt, using the gate, while 
figure l 6 presents the distribution function F(t)) for the same quan- 
tity and practically the same conditions. The interesting fact here 
is the asymmetry of the distribution as first demonstrated by Townsend 
(ref. 3). From the distribution p(t|) the mean values can be computed. 
However, in order to compute the odd moments, that is, du^/dt for 
example, one has to be very sure about the origin of the coordinate T|, 
that is, about the zero position. Since a stationary process is dealt 
with here, evidently 


du 


HpU) dt) 


dt 



P( T]) dT) 


= 0 


Hence one has to obtain first the "center of gravity" of the distribu- 
tion and then compute the moments about this point. For the even moments 
this is not very important and a slight shift of the zero point amounts 
to a negligible error. The details of the computation for the case of 
the 10-channel analyzer are given in the appendix. 

The skewness and flatness factors obtained are shown in the following 
tables . Included also is the factor C defined by 


C 



u 
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for u and du/dt, respectively. Results for d 2 u/dt 2 using the 
10-channel analyzer are also given. (Also see fig. 17.) 

The following tables give a comparison of measurements behind a 
turbulence grid with a mesh size M = 3 /A inch; R^ denotes the Reynolds 

number based on the mesh size, and x/m denotes the dimensionless dis- 
tance from the grid position. 

With gate circuit, R^ = 8 , 100 , and X/M = 62 : 


S f C 


u 

0.002 

2.56 

0.9AA ^«/2 

du/dt 

-0.A2 

3.76 

1.02 jt/2 


With 10-channel analyzer, R^ = 11 , 000 , and X/M = 50 : 


S f C 


u 

0.006 

2.73 

1 r— — 

O.98O y Jt/2 

du/dt 

-O.A 39 

3.80 

0.98^ ^ Jt/2 

d 2 u/dt 2 

0.031 

^•59 

1.0A \j rt/2 


These values agree well with Townsend's measurements. Townsend 
gives for the skewness factor of du/dt, S = - 0 - 39 , hut the measure- 
ments plotted, for example in figure 8 of reference 10, range between 
0.39 and 0.50 at least. The flatness factor of u found here is some- 
what less than Townsend's value of 2.9 (ref. ll) . Hence the present 
measurements do not agree as well with a Gaussian distribution, for 
which f = 3. The same result is borne out by C which should be 

^jt/2 for a Gaussian distribution. Again the distribution of u does 
give a noticeable difference. The symmetry of the u distribution is, 
however, well shown by S, which is zero to a close approximation. 

Joint Probability Density p(|,T|) From Zero Counts 

During the investigation of the possibilities of using the number 

of zeros of u(t) to determine (du/dt) 2 (or A) it was necessary to 
assume that u(t) and du/dt were statistically independent; that is. 
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it was assumed that p(|,T]) = p 0 (|)p 1 (t]). Using the present device with 
the original gate circuit it is possible to check how closely this assump- 
tion agrees with the facts, that is, how accurately (du/dt)^ can be 
measured from the zero counts . 

The number N(£) of | values of u(t), that is, the number of 
times per second that u(t) passes |, is given by (refs. 1 and 12) 



h|p(i,ri) dT i 


(23) 


Denote the probability distribution of u(t) , that is, the average time 
that u(t) spends between £ and £ + d|, as Pq(^) d£ . Evidently 

Pq(|) can be written 


p 0 (i) 



P(l>Tj) dTJ 


Taking the ratio N/p^ one has 


N(S) 

p 0 (i) 


f | Ti|p(6,n) dT] 



p( t > Tl) dT) 


(2k) 


Equation (24) describes the following facts: Pq(^) denotes the 

time t which u(t) spends between £ and £ + d| , and N(|) gives 
the number of passages through £ per unit time; hence the ratio 


P 0 (l) d£ 


T 
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gives the average time of a single passage. Hence the average absolute 
slope of a passage is given by 


In general, |r)|^ will depend upon £ ; however, if p(|,rj) = Pq(£)P;l(t)) 


which is independent of This relation means that, if the probabilities 

of u and du/dt are independent, the absolute mean of du/dt is the 
same at every point of the trace u(t) . 

An example for which this is evidently not true is a simple harmonic 
wave of angular frequency that is 


Here, N(£) = a>/rt is of course independent of However, the time 

spent between £ and | + d£ increases with £ and hence 


then 



(25) 


u(t) = sin a)t 


N(6) 


0 if | — ► 1 



0 


expressing the fact that the mean absolute slope 
at 1=0 and decreases to zero as £ — > 1. 


| T) | has its maximum 
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Without any assumption on the joint probability p(£,T)) the mean 
absolute value of tj can be written 


41 = 


pU, 4) | 4 1 d £ d 4 


p(Ij 4) d £ d 4 


or 


N(t) di 


P 0 (i) d i 


(26) 


From a measured distribution of N(|) and of Pq(S) both the ratio in 

equation (24) and | ff | according to equation (26) can be obtained and 
compared. The result is plotted in figure 18; N/p^ at 1=0 is 
normalized to 1. If u and du/dt were strictly statistically inde- 
pendent N(|)/Pq(£) would be constant and would agree with the ratio 

of the integrals in equation (26). Actually N(i)/p Q (|) is not quite 

constant but shows a definite minimum and maximum, as seen in figure 18, 
and the ratio of the integral’s is less than unity. At large values 
of | the errors in N/pq evidently become very large indeed since the 

ratio approaches 0/0. Still it is evident that N/pq does not vary 
much over the most important range of | near zero where p Q (|) is 
large. The general behavior of N/pq is quite interesting and should 

be studied further. The error in using the value of N(0) to determine 
| tj | is quite small as seen from figure 18 (namely, 1 - O.98 = 0.02). 

There remains one other correction which has to be applied to A 
measurements from zero counts . Evidently the zero counts furnish 

directly only the mean slope of u(t) , that is, | du/dt |. On the other 
hand A is directly related to \| (du/dt) 2 . Hence one has to relate 
| du/dt | to ^ (du/dt) 2. In the simplest formulation one assumes P 1 (’n) 


NACA TN 3037 


25 


to "be Gaussian and then the ratio 


du 

dt 



Now P]_(t]) is not quite Gaussian and hence the relation between 
and ^(du/dt)^ will differ from ^2 /rt . Let 


du/dt | 



Then the formula employed in the measurement of A by counting zeros 
should read 


n 2 (o) = 


u 2 

C 2 2*A 2 


where U denotes the mean speed. For C = fj 2 the formula evidently 
reduces to the one used in reference 1 

N 2 (0) = — (27) 

* 2*2 

The factor C was evaluated from the measured distribution p-^(ti) as 
has been shown in the preceding table: 


I - 


and hence 


n 2 (o) 


0.9^ 


u 

rt 2 A : 


(28) 
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Q 

Consequently in A about 6 percent error results from the use of the 
simple formula, equation (27) • Thus A would come out too large using 
equation (27). This is an error in the right direction since the measure- 
ments of references 1 and 13 consistently give too large values of A as 
calculated from equation (27). The correction due to the weak statistical 
dependency of u and du/dt is, however, in the other direction and the 
combined effect of these two errors is insufficient to explain the dif- 
ferences found in reference 13 between A from zero counts and A 
measured independently. It seems more likely that these differences can 
be traced to small systematic errors in the equipment, for example, insuf- 
ficient resolution. 


Vortex Street 

Figure 19 shows the distribution function, by the 10-channel analyzer, 
obtained in a vortex street behind a cylinder at a Reynolds number of 100. 
The hot-wire output was a very stable, almost triangular wave. A true 
triangular wave would give a straight-line distribution function. The 
rounded peaks of the actual wave form give deviations from a straight 
line near the ends, where the distribution curves inward as in the dis- 
tribution function for the sine wave. This work was done in conjunction 
with the study of wakes (ref. 14) . 


California Institute of Technology, 

Pasadena, Calif., September 29, 1952. 
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APPENDIX 

CALCULATION OF POSITION OF g = 0 AXIS 

Figure 20 is to be used in conjunction with this discussion of the 
calculation of the | = 0 axis. The distance e of the zero axis from 
channel 5 is found by the following method. 

For zero average value of the input signal, 

A 1 + a i = A 2 + a 2 
a 2 — a ^ — A — Ag 

where A^ and Ap are determined from the data by use of Simpson's rule. 

The curve between channels 5 and 6 can usually be accurately approxi- 
mated by a straight line. Therefore, a-|_ and a 2 are computed by using 

the trapezoidal rule as follows: 

a l = |( h l + l l) e 

a 2 = ^(^2 + ^ 2 )^ “ € ) 

2 ^2 - 0*2^ = ^2 *” (^1 + ^2)^ 

= 1^2 ^2 “ ^h-2. ^ ^~) ^ 

But, 

(l - hg - hje 


= 1-1! 


*1 = h i + 
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Substituting, 


2^a 2 - a-^ = hg + 1 - + (b.-^ + hg - l)e - ^h-]_ + h 2 ^e - € 

2 ^2. “ Ag^ = hg - h^_ + 1 - 2e 

Therefore , 


hg - + 1 

€ = 

2 
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(b) Square-wave signal from gate circuit. 



(c) Signal from pulse generator. 



and pulse-generator circuits. 
Figure 1.- Principle of gate circuit. 
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(a) Fluctuating voltage. 



(b) Unmodulated pulses. 




1 

(d) Wave pattern from discriminator. 

Figure 2.- Principle of 10-channel analyzer. 
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REGION OF INTEGRATION 
IN EQUATION (IS) 


Figure 3*- Typical distribution function. 
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© POTENTIOMETER FOR SETTING \ VALUES 
© TRIGGERING ADJUSTMENT FOR CIRCUIT (l) 

© TRIGGERING ADJUSTMENT FOR CIRCUIT (Z'i 
® OUTPUT EQUALIZATION CONTROL 

© METER TO ENABLE CURRENT IN ©TO BE KEPT CONSTANT 

Figure 4.- Double -trigger circuit of gate. (K indicates thousand; 

WW , wire-wound resistor.) 
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Figure 5 .- Circuit diagrams of inverter and mixer. 


(K indicates thousand.) 



Figure 6.- Block diagram of apparatus. 
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[— 500 p-SEC -j 


L-802I4.5 

Figure 7*- View of a single pulse. 
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2,000 PULSES/ SEC 


L-802IJ.6 

Figure 8.- Differential amplifier output for a sine-wave input. 


38 


NACA TN 3037 



1 

1 1 

1 i 

. 1 1 

1 1 

1 1 1 

1 1 1 

1 1 1 

1 1 1 

till 

1 1 1 

1 1 

III 

1 1 1 

1 1 1 

.1111 

1 1 1 

1 1 I 

Hill 

1 1 1 

1 1 1 

1 1 1 1 1 1 

1 1 


2,000 PULSES /SEC 


L-802J4.7 

Figure 9*- Modulator output for a sine-wave input. 


♦ 300 V 



T 7 - fl A H 6 
T a - 8AL5 


Figure 10.- Circuit diagram of pulse amplitude modulator. 
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Figure 11.- Pulse height discriminator. 
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Figure 13*- Distribution of u in isotropic turbulence (corrected for gap 
width). R m = 8,100, X/M = 62; C = 0.944/«/2, S = 0.002, f = 2.56. 
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Figure Ik . - Distribution function for u(t) in isotropic tu rbule nce . 
% = 11,000, X/M = 50, M = 1.68 centimeters; C = 0.980\ZrT /2, 

S = 0.006, f = 2.73. 
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Figure 15.- Distribution of du/dt in isotropic turbulence. R M 
X/M = 62; C = 1.02\/n/2, S = -0.42, f = 3 . 76 . 
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Figure 16 .- Distribution function for du/dt in isotropic turbulence. 
= 11,000, X/M = 50, M = 1.68 centimeters; C = 0.984/n /2, 

s = -0A39, f = 3.80. 
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Figure 17 •- Distribution function for d 2 u/dt 2 in isotropic tur bulence. 

% = 11,000, X/M = 50, M = 1.68 centimeters; C = 1.04j/jt/2, S = 0.031, 

f = 4.59. 
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Figure 18.- Plot of NCO/P-lCi). 



* 

Figure 19*- Distribution function for u(t) in vortex street behind 
circular cylinder of diameter d at a Reynolds number R = 100. 
x/d = 6, Y/d = 1.3; C = 1.15, S = 0.08, f = 1.7. 
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Figure 20.- Calculation of | = 0 axis. A^_, A 2 , a^_, a 2 > areas; 

Z2, h-^, h2, €, lengths. 
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